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ABSTRACT 



We present unique spectroscopic radio observations of comets C/2002 X5 (Kudo- Fujikawa), C/2002 VI (NEAT), and 
C/2006 PI (McNaught), which came within r h « 0.2 AU of the Sun in 2003 and 2007. The molecules OH, HCN, 
HNC, CS, and CH3OH were detected in each of these comets when they were exposed to strong heating from the 
Sun. Both HC 3 N and HCO+ were detected in comets C/2002 X5 and C/2006 PI, respectively. We show that in these 
very productive comets close to the Sun screening of the photodissociation by the Sun UV radiation plays a non- 
negligible role. Acceleration of the gas expansion velocity and day-night asymmetry is also measured and modeled. 



The CS photodissociation lifetime was constrained to be about 2.5 x 10 



at rn = 1 AU. The relative abundances 



are compared to values determined from more distant observations of C/2002 X5 or other comets. A high HNC/HCN 
production-rate ratio, in the range 10-30% between 0.5 and 0.1 AU from the Sun, is measured. The trend for a significant 
enrichment in CS in cometary comae (CS/HCNoc r^ ' 8 ) is confirmed in all three comets. The CH3OH/HCN production 
rate ratio decreases at low ru- The HC3N/HCN production rate ratio in comet C/2002 X5 is four times higher than 
measured in any other comet. 

Key words. Comets: general - Comets: individual: C/2002 X5 (Kudo- Fujikawa), C/2002 VI (NEAT), C/2006 PI 
(McNaught) - Radio lines: solar system - Submillimeter 



1. Introduction 

The composition of cometary nuclei is of strong interest 
in understanding their origin. Having spent most of their 
time in a very cold environment, these objects should not 
have evolved much since their formation. Thus, their com- 
position provides clues to the composition in the outer 
regions of the Solar Nebula where they formed. The 
last two decades have proven the efficiency of microwave 
spectroscopy in investigating the chemical composition 
of cometary atmospheres. About 20 different cometary 
molecules have now been iden tified at radio wavelengths 
(jBockelee-Morvan et aill2004ari . 

In this paper, we extend our investigations of the com- 
position of cometary atmosphe r es from rad i o observations 
(iBockelee-Morvan et ail l2004eJ: iBiver et all l2002al l2006aL 
|2007y) to three comets observed in 2003 and 2007. These 
observations were designed to measure the molecular abun- 
dances of comets approaching close to the Sun to investi- 
gate how the strength of the solar heating of both the comet 
nucleus and of its environment affects the coma composi- 
tion. Previous observations have suggested that the rela- 
tive production rat es of severa l molec ules vary with he- 
liocentric distance ([Biver et all 1200 6a). The passages of 
comets C/2002 X5 (Kudo-Fujikawa), C/2002 VI (NEAT), 
and C/2006 PI (McNaught) provided us with the opportu- 
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nity to measure molecular abundances at heliocentric dis- 
tances (r/j) between 0.1 and 0.25 AU, about one order of 
magnitude smaller than usual. This study is complementary 
to the long-term monitori ng of comet C/1995 01 (Hale- 
Bopp) (jBiver et all l2002bl) , which provided information on 
the outgassing of a comet between 0.9 and 14 AU. 



Opportunities to plan observations of comets passing 
within 0.2 AU from the Sun are rare. For example, comet 
C/1998 Jl (SOHO) was discovered too late to establish 
an accurate ephemeris around perihelion time. In addi- 
tion, Sun-grazing comets often do not survive and even 
disintegrate before reaching perihelion, making observing 
plans extremely difficult. The last opportunity was comet 
C/1975 VI (West). The observations reported here are 
unique, and required to develop specific observing strate- 
gies and analyses. 



The organization of the paper is as follows. In Sect. 
2, we present the observations of comets C/2002 X5 
(Kudo-Fujikawa), C/2002 VI (NEAT), and C/2006 PI 
(McNaught) performed with the 30-m telescope of the 
Institut de Radioastronomie Millimetrique (IRAM) and the 
Nancay radio telescope. In Sects. 3 and 4, the analysis of 
these observations is presented. A summary follows in Sect. 
5. 
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2. Observations 

Owing to the late discovery or assessment of their activity, 
the three comets were observed as targets of opportunity 
using the IRAM 30-m and Nancay radio telescopes. This 
was possible because these telescopes do not have tight solar 
elongation constraints. The small solar elongation (< 10° at 
r/i < 0.2 AU) resulted in limited observing support from 
other observatories. 



2.1. C/2002 X5 (Kudo-Fujikawa) 

Comet C/2002 X5 (Kudo- Fujikawa) was discovered visu- 
ally at m v — 9 on 13-14 December 2002 by two Japanese 
amat eur astronomers, T. Kudo and S. Fujikawa (jNakanol . 
12002ft . At about 1.2 AU from the Earth and the Sun at 
that time, it was then a moderately active comet. It passed 
perihelion on 29 January 2003 at a perihelion distance q = 
0.190 AU from the Sun. 

One of the main goals of the observations of comet 
C/2002 X5 at the IRAM 30-m was to measure the evo- 
lution of the HNC/HCN production rate ratio as it ap- 
proached the Sun. A first observing slot was scheduled on 4- 
5 January 2003, at = 0.8 AU, but the weather prevented 
any observation. Observations performed on 13 January 
(r/j = 0.55 AU) only partly succeeded because of techni- 
cal problems. Most data were acquired on 26.5 January, 2 
days before perihelion, at r/, = 0.21 AU. At that time, the 
comet was only observable visually by the C3 coronagraph 
aboard the SOlar Heliosph eric Observatory (SOHO) — the 
solar elongation was 5.5° (|Bout et all . 120031: iPovich et all . 
2003) . The pointing of the comet was a real challenge, since 
the ephemeris uncertainty was expected to be on the order 
of the beam size (10-20"). Since at = 0.2 AU molecu- 
lar lifetimes are typically shorter than an hour, hence pho- 
todissociation scale lengths are smaller than the beam size, 
accurate pointing was required. From coarse mapping, we 
found the comet about 10" south of its predicted position 
(Fig. [U Table [[]). The last observing run at IRAM took 
place on 12 March 2003 (r h = 1.2 AU, Fig. [3]), when the 
comet was receding from the Sun and about 100 times less 
productive. A log of the observations and measured line 
areas are given in Table [2] Sample spectra are shown in 
FigsHU 

To monitor the water production rate, observations of 
OH at 18-cm using the Nangay radio telescope were sched- 
uled on a daily basis from 1 January to 10 April 2003. The 
OH lines were only detected when the comet was between 
1 and 0.4 AU from the Sun, inbound and outbound. 



2.2. C/2002 VI (NEAT) 

Comet C/2002 VI (NEAT) was discovered on 6 November 
2002 by the Near Earth Asteroid Tracking (NEAT) pro- 
gram telescop e on the Haleakala summit of Maui, Hawaii 
(|PravdoLl2002D . At that time, it was a relatively faint object 
(m v = 17). Comet C/2002 VI (NEAT) brightened rapidly 
after its discovery, becoming a potentially interesting tar- 
get for observations at perihelion on 18 February 2003 at q 
= 0.099 AU. The orbital period was estimated to be 9000 
years (Nakano note NK965Q), suggesting that it has sur- 
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Fig. 1. Coarse map of the HCN J = 3-2 emission from 
comet C/2002 X5 (Kudo- Fujikawa) observed with the 
IRAM 30-m. Positional offsets in arcsec are with respect 
to the final ephemeris (Table Q]). All four spectra are plot- 
ted at the same intensity and velocity scales, provided in 
the upper right box. The projected direction of the Sun is 
also given. 
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Fig. 2. CH 3 OH lines at 157 GHz observed with the IRAM 
30-m telescope in comet C/2002 X5 (Kudo- Fujikawa) on 
26.5 January 2003. The vertical scale is main-beam bright- 
ness temperature. The horizontal scales are the rest fre- 
quency (upper scale) or the Doppler velocity in the comet 
rest-frame, relative to the (6o — 6_i)E line frequency (bot- 
tom scale). 



vived a close passage to the Sun at its previous perihelion 
and might be expected to do so again. 

The observations at IRAM were undertaken on 16 and 
17 February 2003 (r h = 0.13-0.11 AU), i.e., just one day be- 
fore perihelion. As for comet C/2002 X5 (Kudo- Fujikawa), 
the observations were challenging because of the lack of sup- 
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Table 1. Observing conditions and ephemeris errors 



Comet 


UT date 


Th. 


A 


Phase 


O - 


/-i 1 


Cf 






[yyyy/ mm / dd - d J 


r A T Tl 

[AUJ 


r A T Tl 

[AUJ 


[ 1 




ri~\ _ _ 

ouec. 


rri a 

oKA 


ri~\ _ _ 

ouec. 


C/2002 X5 


2003/01/13.6 


0.553 


1.032 


69.3 






+1.6" 


-2.2" 




2003/01/26.5 


0.214 


1.172 


26.1 


1.5" 


-11" 


+0.9" 


-10.3" 




2003/03/12.7 


1.184 


1.096 


50.5 






+3.0" 


+0.9" 


C/2002 VI 


2003/02/16.6 


0.135 


0.978 


90.3 


-1.9" 


-3.5" 


-3.9" 


-6.6" 




2003/02/17.6 


0.107 


0.986 


87.7 


-1.5" 


-4.5" 


-2.5" 


-8.8" 


C/2006 PI 


2007/01/15.6 


0.207 


0.817 


140.0 


+1" 


+30" 


-1.2" 


+27.1" 




2007/01/16.6 


0.229 


0.822 


129.3 


+1" 


+28" 


+0.6" 


+25.7" 




2007/01/17.6 


0.256 


0.833 


119.1 


+0" 


+23" 


+2.0" 


+23.5" 



1 O, C u , and Cf denote the position of the HCN peak emission determined from mapping, the used computed ephemeris, and 
the final (reference) ephemeris, respectively. Orbital elements are given for C/2002 X5 (MPEC 2003- A41, 2003- A84, and 2003-D20 
used, JPL#40 2 reference), C/2002 VI (MPEC 2003-C42 used, JPL#34 reference), C/2006 PI (JPL#15, JPL#25 reference). The 
uncertainty in reference orbit is 2". 



JPL HORIZONS ephemerides: http://ssd.jpl.nasa.gov/7ephemerides 
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Fig. 3. HCN J = 3-2 line observed with the IRAM 30- 
m telescope in comet C/2002 X5 (Kudo-Fujikawa) on 12.7 
March 2003. The vertical scale is the main beam brightness 
temperature. The horizontal scale is the Doppler velocity 
in the comet rest-frame. 
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porting optical astrometry, the comet being 8-6° away from 
the Sun and only seen by the SOHO coronagraph when it 
became as bright as magnitudes —1 to —2. A coarse map of 
the HCN J=3-2 line at IRAM revealed the comet to be at 
about 5.5" (half a beam, Fig. [51 Table[JJ) from its expected 
position. Observational data are given in Table [31 

OH 18-cm observations were obtained nearly every day 
from 31 December 2002 to 25 March 2003. The comet 
was only detected at > 0.4 AU from the Sun, though 
the water production rate (Qh 2 o) likely exceeded 2 x 10 30 
molec. s _1 at perihelion (Sect.3). 

2.3. C/2006 PI (McNaught) 

C/2006 PI was discover ed by Robert Mc Naught at m v = 
17.3 on 7 August 2006 (|McNauehtJ . l2006h as it was at 3.1 
AU from the Sun. The geometry was very unfavorable for 
observing this comet as it approached the Sun. The comet 
was basically lost in the glare of the Sun in November and 
December 2006 from r h = 1.5 to 0.5 AU. At r h > 0.5 AU, 
its intrinsic brightness was similar to that of comet C/1996 



Fig. 4. HC 3 N, CS, HCN, CH 3 OH (sum of the 3 brightest 
lines at 157 GHz) and HNC lines observed with the IRAM 
30-m telescope in comet C/2002 X5 (Kudo-Fujikawa) on 
26.5 January 2003. The spectra are plotted with the same 
horizontal scale to compare the line widths. The intensities 
of the spectra of HCN and HNC were divided by 12.5 and 
2.5, respectively. Spectra are shifted vertically for clearer 
viewing. 



B2 (Hyakutake). It brightened rapidly in early January 
2007 to peak at m v ~ -5 and becam e visible to the naked 
eye in broad daylight (jGreenl . 120071 ). It passed perihelion 
on 12 January 2007 at q = 0.17 AU. Comet C/2006 PI 
(McNaught) was the brightest and most productive comet 
sin ce C/1965 SI (Ike ya-Seki). The absence of an ion tail 
led lFulle et al.l (|2007|) to argue that, because of a very high 
outgassing rate, the diamagnetic cavity was so large that 
ions were photodissociated before reaching the region where 
they could interact with the solar wind. A week after peri- 
helion, C/2006 PI displayed a fantastic dust tail with many 
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Fig. 5. Scales as in Fig. [TJ coarse map of the HCN J = 
3-2 emission in comet C/2002 VI (NEAT) with the IRAM 
30-m on 17.5 February 2003. But in this case, spectra are 
plotted at positions relative to the ephemerides used during 
the observations. 
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Fig. 6. CS, HCN, CH 3 OH (sum of the 8 brightest lines at 
252 GHz), and HNC lines observed with the IRAM 30- 
m telescope in comet C/2002 VI (NEAT) on 16.5-17.6 
February 2003. The spectra are plotted with the same hor- 
izontal scale to compare the line widths. The intensities of 
the spectra of HCN and HNC were divided by 10 and 2, 
respectively. Spectra are shifted vertically for clearer view- 
ing. 
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Fig. 7. On-the-fly map of the HCN J = 3-2 emission in 
comet C/2006 PI (McNaught) obtained with the IRAM 
30-m on 15.65 January 2007 (4 min. integration time). The 
de-pixelized image shows the line area between —3 and +3 
km s , contour levels are 2 and 4er rms . The cross marks the 
comet position as computed from the reference ephemeris. 
Before gridding to beam/2 sampling, a linear spatial base- 
line was subtracted from each on-the-fly subscan to remove 
total power fluctuations of mainly atmospheric orgin. 



striae curving around 1/3 of sky at a mean distance of 
20° from the Sun. 

The IRAM observations were performed on 15, 16, and 
17 January 2007 (r h = 0.21-0.25 AU). The comet ephemeris 
was expected to be possibly wrong by up to 1' (correspond- 
ing to the 3-g uncertainty fro m the JPL's HORIZONS 
ephemeris (jGiorgini et all [1996)). An on-the-fly map of the 
HCN J(3-2) line on 15.6 January UT with the IRAM 30- 
m HERA array of receivers showed that the comet was 
30" north of its predicted position (Fig. [7]). Orbit updates 
later confirmed the observed offset (Table [T]). Most observa- 
tions consisted of five-point integrations at and 6" offsets 
in RA and Dec to pinpoint the maximum emission (e.g., 
Fig.[S]). On 17.5 January, the observations were affected by 
strong anomalous refraction due to the low elevation of the 
comet (21-17°). The effect was estimated to be equivalent 
to a mean pointing offset of up to 12", implying a loss of 
90% of the signal. Observational data are summarized in 
Table H 

OH 18-cm observations were performed daily between 
10 and 20 January. The comet was detected intermittently. 
OH 18-cm lines had never been detected that close to the 
Sun (0.17 AU) in a comet before. 



3. Analysis of Nancay data - HoO production rates 

The characteristics of the Nangay radio tel escope and the 
OH ob servations of comets may be found in lCrovisier et all 
(120021) . The OH 18 -cm lines usually observed in comets are 
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Fig. 8. HDO line at 225.896 GHz observed with the IRAM 
30-m telescope in comet C/2006 PI on 16.6 January 2007. 
The vertical scale is the main-beam brightness temperature. 
The horizontal scale is the Doppler velocity in the comet 
rest-frame. 
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Fig. 9. Coarse map of the HCN J(3-2) emission from comet 
C/2006 PI (McNaught) on 16.5 January 2007, displayed 
as in Figs [T] and O with positions relative to the final 
ephemerides (Table [T]). 



maser lines pumped by UV solar radiation ()Despois et all 
Il98lh . However, for high water-production rates, the maser 
emission is quenched by collisions in a large part of the 
coma, and most of the signal may originate from ther- 
mal emission. The radius of the region where the maser 
emission is quenched is estimated to be 50000^-^029 km, 
where Q29 is the water product ion rate in units of 
10 29 molec. s" 1 (|Oerard et allll998fl . For production rates 
above 6 x 10 29 molec. s _1 and the heliocentric distances 
where the comets were observed (r^ < 0.3 AU), most OH 
radicals dissociate within the quenching zone. Hence, the 
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Fig. 10. CS, HCN, CH3OH (sum of the 4 brightest lines at 
157 GHz), and HNC lines observed with the IRAM 30-m 
telescope in comet C/2006 PI (McNaught) on 16.5 or 17.6 
January 2007. Spectra have been plotted at the same hori- 
zontal scale to compare the line widths. The intensity scale 
is reduced by 10 and 2 for HCN and HNC, respectively, and 
shifted vertically for better viewing. 



signal is dominated by OH thermal emission. The OH pro- 
duction rates or upper limits are given in Tableland sam- 
ple spectra are provided in Fig. [3] 

Table [5] also includes water-production rate measure- 
ments from other instruments for comet C/2002 X5. The 
H2O line at 557 GHz was observed in th is comet with 
Odin during March 2003 (jBiver et alll2007ah . Observations 
were conducted with the SOHO coronagraph spectrom- 
eter on the perihelion date: from H Lyman a measure- 
ments, the peak outgas sing rate of water is estimated to 
3 x 10 30 molec. s _1 (ICombi et all l2008ft . According to 



iPovich et all (|2003ft and lBout et al 



20011 7 this comet dis- 



played a strongly variable activity around perihelion with 
a ^two day period. During the perihelion period, the 
comet was not or only very marginally detected at Nangay. 
Averaging the data obtained during the five days around 
perihelion, we obtain a possible 4-cr detection (subject to 
baseline uncertainties) suggesting a production rate of ~ 
6x 10 30 molec. s" 1 (Table©. This is higher than the SOHO 
estimate, but still within the same order of magnitude. 
Power laws fitted to Nangay and Odin data (Fig. [T5| yield 
values in the same range 1.5-5.4x 10 30 molec. s _1 at 0.214 
AU from the Sun. A mean value of 3.5xl0 30 molec. s _1 is 
adopted. 

For C/2002 VI (NEAT), the pre-perihelion observations 



can be fitted by a power law (0.33±0.03x 10 2 



-1.95±0.23 



), 



which extrapolates to Qh 2 o = 1-6 x 10 30 and 2.5 x 10 30 
molec. s _1 for 16 and 17 January 2003, respectively, con- 
sistent with the upper limit determined for these days. 
Extrapolation of SOHO-SWAN measurements yield values 
that are four times higher but indicative of abundances rel- 
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ative to water that are abnormaly low for all molecules. We 
assume that <5h 2 o = 2.0 and 2.5 x 10 30 molcc. s _1 , respec- 
tively, these values being more compatible with the day- 
to-day variations observed in millimeter spectra of other 
molecules. 

As for C/2006 PI (McNaught), the water production 
rates for 13 and 19 January deduced from the Nangay 
data agree with the estimated HDO production rate from 
the marginally detected line at 225.9 GHz (Fig. [5]) and 
the HDO/H20= 6 x 10~ 4 ratio measured in comets. The 
total production rate varied from w 40 x 10 30 down to 
ps 5 x 10 30 molec. s _1 at that time (Sect. 3). 
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The estimated visual magnitudes of these comets at 
perihelion were mi ss +3, mi —1.5, and mi 
-5 for C/2002 X5, C/2002 VI, and C/2006 PI, respec- 
tively. Correcting for the —2 magnitude surge i n bright- 
ness of C/2006 PI caused by forward scattering (|Marcusl . 
l200l . and using the correlation law between heliocen- 
tric magnitude and Qh 2 o of iJorda et~ai1 (|2008l ). this 
would imply water production rates of 1, 11, and 26xl0 30 
molec. s _1 respectively. The comparison to measured pro- 
duction rates in Table [5] suggests that C/2002 VI was 
a more dusty comet than the two others since, unlike 
C/2002 X5 and C/2006 PI, the actual outgassing rate being 
much lower than the value (11 x 10 30 molec. s _1 ) inferred 
from visual magnitudes. 



4. Analysis of IRAM data 

Molecular production rates were derived usi ng models of 
mole c ular excitati on and radiation transfer (jBiver et all 
119991 120001 l2006al) . The excitation model of CH 3 OH was 
updated. The computation of the partition function con- 
siders now the first torsional state, which is populated sig- 
nificantly in the hot atmospheres of the comets studied in 
this paper. For these productive comets observed at small 
Th , a significant fraction of the observed molecules are pho- 
todissociated before leaving the collision-dominated region. 
Thus, the determination of the gas kinetic temperature 
(which controls the rotational level populations in the col- 
lision zone) and strong constraints on molecular lifetimes 
were essential. As shown below, line shapes and brightness 
distributions obtained from coarse maps provide informa- 
tion on the gas outflow velocity and molecular lifetimes. 

4.1. Gas temperature 

Table [B] summarizes rotational temperatures deduced from 
relative line intensities. Using our excitation models, we 
then constrained the gas tem perature (Ta ble [5]) following 
the methods outlined in, e.g., iBiver et al.l (|1999|) . Several 
data are indicative of relatively high temperatures, which 
are difficult to measure because the rotational population 
is spread over many levels making individual lines weaker. 
Hence, the uncertainties in derived values are high. We note 
that the values derived for C/2002 X5 from CH 3 OH lines 
do not reflect the marginal detection of some lines. The in- 
consistency between the different measurements, especially 
from CH 3 OH and HCN lines in comet C/2006 PI, is pos- 
sibly related to differences in collision cross-sections and 
temperature variations in the coma. 

Gas temperature laws as a function were ob- 
tained for comets C/1995 01 ( Hale-Bopp 1 ) (IBiver et all 
I2002bh . C/1996 B2 (Hyakutake) (IBiver et all Il999h. and 



153P/Ikeya-Zhang (jBiver et alll2006al) . The laws (100< 
for Hale-Bopp and 60r^ ' 9 for other comets on average), 
extrapolate to T = 200-1000 K for r h = 0.25-0.11 AU. 
Photolytic heating increases with decreasing and the in- 
creasing production rate of water. Hence, we expect higher 
temperatures for the mo st productive comet C/2006 PI. 
ICombi and Smyth! (|1988ft predicted a maximum temper- 



ature on the order of 500-900 K for comet Kohoutek at 
0.25-0.14 AU, but temperatures are generally below 300 K 
at distances from the nucleus where molecules are not pho- 
todissociated. Given also that the HCN J(3-2) V2 = 1 line 



at 265852.709 MHz is not detected, we estimate that tem- 
peratures are lower than 300 K below r/j=0.25 AU. 

For C/2002 X5, we adopted the gas kinetic tempera- 
ture values of 50 K, 180 K, and 40 K for mid- January, end 
of January, and mid-March all in 2003, respectively. For 
C/2002 VI and C/2006 PI, we used 150 K and 300 K, re- 
spectively. We later discuss (in Sect. 5.1) the influence of 
the assumed temperature on the inferred production rates. 

4.2. Screening of photolysis by water molecules 

Given high water-production rates (> 10 30 -10 31 
molec. s^ 1 ), self-shielding against photodissociation 
by solar UV radiation is significant for these comets. As 
a consequence, molecules can have longer lifetimes on the 
night side, and reduced H2O photolysis may limit gas 
acceleration. We studied the screening of photolysis using 
simplified assumptions: isotropic outgassing at a constant 
expansion velocity and an infinite lifetime for the screening 
molecules (i.e., we assumed that their scale length is larger 
than the s iz e of t he optically thick r egion). Accor d ing to 
iNee fc Led (Il985h . when considering iLee fc Sutol (|l986h 
and iLed (|l984h ~ the OH photoabsorption cross-section 
does not differ much from the H2O one: it peaks close 
to the Lyman a wavelength and is slightly stronger (1-2 
times). On the other hand, the cometary hydrogen may 
not absorb significantly the Solar Lyman a spectral line. 
The comet Lyman a line is too narrow (velocity dispersion 
of 20 km s -1 ) to absorb significantly the solar emission 
(«180 km s _1 line width). Consequently, if the sum of the 
scalclcngths of H2O and OH is larger than the distances 
from the nucleus considered hereafter, the infinite lifetime 
assumption should not underestimate the screening effect. 

In cylindrical coordinates, with the vertical z-axis point- 
ing towards the Sun, a point in the coma has coordinates (9, 
p = rsin(</>), z r = rcos(</>)), where r is the distance to the 
nucleus and <f> the co-latitude angle. The photodissociation 
rate of a molecule M, characterized by its photodissociation 
absorption cross-section ctm(A), in a radiation field F (A) (in 



photons m 



nm ) is given by 



(1) 



The problem is symmetric around the z-axis. At (p, z r ), the 
solar flux at wavelength A will be attenuated to 



F(\,p,<j>) = F (A)exp[-r(A,p»] 



(2) 



because of the optical thickness along the comet-sun axis. 
We consider water as the major molecule responsible for the 
opacity, which is connected to its absorption cross-section 
ch 2 o(A) by 



t(X,P, <p) = o-H 2 o(A) / nn 2 o(p,z)dz 



(3) 



The photons absorbed are mostly those responsible for pho- 
todissociation and not fluorescence. Photodissociation of 
molecules takes place for A < 200 nm and a significant 
solar UV field (A > 80 nm). Using the assumptions for the 
density, we can integrate the density along the z-axis to 
obtain 
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Table 5. Water production rates 



UT date 


< r h > 


Observatory 


Line intensity 


Qh 2 o 


Ref. 


[mm/dd.dd] 


[AU] 


and line 


[mjy km s _1 ] 


[molec. s -1 ] 




C/2002 X5 (Kudo- Fujikawa): (2003) 


01/12.5-14.5 


0.55 


Nangay OH 18-cm 


-156 ±7 


1.5 ± 0.1 x 10 2y 




01/15.5-19.5 


0.45 


Nangay OH 18-cm 


-91 ±5 


1.8 ± 0.1 x 10 29 




01/26.50 


0.214 


Nangay OH 18-cm 


47 ± 15 


85 ± 27 x 10 29 




01/25.5-30.5 


0.205 


Nangay OH 18-cm 


« 32 ±8 


w 64 ± 16 x 10 29 




01-02 


< 1 


SOHO-SWAN 




0.56 x 10 29 x r^ 2 


[3] 


01/27.92 


0.195 


SOHO-UVCS Lyman a 




5.5 x 10 29 


[1] 


01/28.13 


0.194 


SOHO-UVCS Lyman a 




9.0 x 10 29 


[1] 


01/28.79 


0.190 


SOHO-UVCS Lyman a 




13.0 x 10 29 


[1] 


01/29.17 


0.190 


SOHO-UVCS Lyman a 




27.0 x 10 29 


[1] 


03/12.5 


1.178 


Odin H 2 557GHz 


879 ± 79 1 


0.21 ±0.01 x 10 29 


[2] 


C/2002 VI (NEAT): (2003) 


01/00-07 


1.30 


Nangay OH 18-cm 


-25 ±4 


0.15 ±0.02 x 10 29 




01/08-14 


1.15 


Nangay OH 18-cm 


-45 ±4 


0.24 ± 0.02 x 10 29 




01/15-19 


1.01 


Nangay OH 18-cm 


-81 ±5 


0.46 ± 0.03 x 10 29 




01/21-25 


0.87 


Nangay OH 18-cm 


-110 ± 6 


0.55 ±0.03 x 10 29 




01/26-30 


0.73 


Nangay OH 18-cm 


-93 ±5 


0.59 ±0.03 x 10 29 




02/01-05 


0.59 


Nangay OH 18-cm 


-85 ±6 


0.70 ±0.06 x 10 29 




02/06-09 


0.45 


Nangay OH 18-cm 


-66 ±7 


1.7 ±0.2 x 10 29 




02/16.50 


0.136 


Nangay OH 18-cm 


-2± 13 


< 150 x 10 29 




01-03 


< 1 


SOHO-SWAN 




0.94 x 10 29 x r~ 21 


[3] 


C/2006 PI (McNaught): 


(2007) 








01/12.51 


0.171 


Nangay OH 18-cm 


45 ±23 


< 270 x 10 29 




01/13.51 


0.173 


Nangay OH 18-cm 


111 ± 15 


400 ± 60 x 10 29 




01/17.52 


0.254 


Nangay OH 18-cm 


-9 ±23 


< 100 x 10 29 




01/19.52 


0.311 


Nangay OH 18-cm 


128 ± 26 


80 ± 20 x 10 29 





1 : l ine integrated inten sity i n mK km s 1 : 

[l]: lPovich et~"aH f200l ); [2]: lBiver et all (|2007al ): [3]: ICombi et~aTl ((20081 ) 

Table 6. Gas temperature measurements 



UT date 
[mm/dd.dd] 


<r h > 
[AU] 


offset 

["] 


Lines 


Rotational temperature 
[K] 


Gas temperature 
[K] 


C/2002 X5 (Kudo- Fujikawa): (2003) 


01/13.61 


0.553 


3.5 


CH 3 OH 157 GHz 


38 ±10 


30 ±10 


01/26.54 


0.213 


3.0 


CH3OH 157 GHz 


320l- 


1 7C+00 


01/26.54 


0.213 


1.6 


CH 3 CN 147 GHz 


1500±~ 05 


> 190 


01/26.54 
03/12.72 


0.213 
1.184 


1.0 
3.0 


HCNJ(3-2),w 2 = l/v 2 = 
CH 3 OH 157 GHz 


< 245 
851" 


< 530 1 

,u -20 


C/2002 VI (NEAT): (2003) 


02/17.0 
02/17.25 
02/17.3 


0.120 
0.116 
0.110 


1.9 
1.6 
1.6 


HCNJ(3-2),« 2 = l/v 2 = 
CH3OH 252 GHz 
CH3OH 252 GHz±20 3 A + ~ 


< 285 

83^27 

228 +92 


< 574 1 
220t 99 


02/17.6 


0.108 


1.8 


HCNJ(3-2)/(l-0) 




80111 


C/2006 PI (McNaught): 


(2007) 








01/16.56 


0.228 


3 


HCNJ(3-2)/(l-0) 


^»»_204 


> 375 2 


01/16.56 
01/16.55 


0.228 
0.229 


5 

5.9 


HCNJ(3-2)/(l-0) 
CH3OH 157 GHz 


> 168 2 


> 600 2 
1 10+°° 


01/16.56 


0.229 


4.9 


CH3OH 157 GHz 


88l? 2 


7f) +20 
,u -10 


01/17.5 
01/17.55 


0.255 
0.255 


6.5 
6.7 


HCNJ(3-2),w 2 = l/v 2 = 
CH3OH 157 GHz 


< 278 
> 130 


< 677 1 
> 90 


01/17.55 


0.255 


6.9 


CH 3 CN 147 GHz 


78l 2 ° 8 6 


80± 3 i° 



1 A significant fraction of the molecules are outside the collisional region and a higher gas temperature is needed to populate the 
v — 1 level inside the collision zone. 

2 A 10% calibration uncertainty in each line is assumed. 
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C/2002 X5 (Kudo-Fujikawa) 
030112-030114 




C/2002 X5 (Kudo-Fujikawa) 
030115-030119 





C/2002 X5 (Kudo-Fujikawa) 
03/01/26.50 



C/2002 X5 (Kudo-Fujikawa) 
030124-030128 



Velocity [km s '] 



C/2002 V1 (NEAT) 
021231-030107 

Mr 




-5 


5 


C/2002 VI (NEAT) 
■ 030126-030130 


-5 


5 


C/2006 P1 (McNaught) 
07/01/12.51 




ft 



C/2002 V1 (NEAT) 
030108-030114 




C/2002 V1 (NEAT) 
030201 -030205 



C/2006 P1 (McNaught) 
07/01/13.51 



Velocity [km s '] 



C/2002 V1 (NEAT) 
030115-030119 





C/2002 VI (NEAT) 
030206-030209 




C/2006 PI (McNaught) 
07/01/17.52 



Velocity [km s '] 



C/2002 VI (NEAT) 
030121-030125 





C/2002 V1 (NEAT) 
03/02/1 6.50 




C/2006 P1 (McNaught) 
07/01/19.52 



Velocity [km s ] 





Fig. 11. Nancay spectra of OH 18-cm lines (average of 1665 and 1667 MHz transitions) corresponding to observations 
listed in Table [SI 



n K2 o(p 1 z)dz 



Qh 2 



z r —p/ tan(0) 

47T7J cxp (p 2 + Z 2 ) 

<9h 2 o <f> 



dz 



iTTV exp r sin(</)) 



Hence at any point in the coma, we estimate the opacity 

40 



T(\,r,4>) =< r(A,r) > 



7T 2 sin(0) 



r sin(0) ' 



(4) 



where the average value of the opacity over 4ir steradians 
at the distance r from the nucleus is 



<r(A,r) >=C{\)~ 
4 r 



with 



C(A) = 



0H 2 o(A)Qh 2 



47TW, 



(5) 



(6) 



cxp 



The surface corresponding to an opacity t(A, r, 4>) — 1-0 is 
defined by 



rthick{(f) = C(A) 



o"h 2 o(A)Qh 2 o 



sin((j)) 



Attv, 



cxp 



sin(0) 



(7) 



According to Eq. [3 the size of the optically thick region 
tends to infinity in the anti-sunward direction (4> = tt). 
It has a finite length determined by the apparent size of 
the Sun (2.7° at r h = 0.2 AU). This region is plotted in 
Fig. HU for the three comets, considering only absorption of 
Lyman a photons by water m olecules with a cro ss-section 
an 2 o(Ly a) = 15 x 1CT 22 m 2 (jLewis et all Il983l) . 

We next consider HCN photodissociation, but similar 
equations can be established for other molecules. The ef- 
fective HCN photodissociation rate at a point (r,cf>) in the 
coma can then be derived from Eq. [1] 



C(A)0 



(8) 



/?hcn {f, — / cTHCN(A)F (A)exp . . 

where <7hcn(A) is the photodissociation cross-section for 
HCN. The integration can be divided over several wave- 
length intervals, corresponding to the different absorption 
bands of HCN and H 2 



A; 



/?HCN,Ai (r, 4>) 



= Axhcn x ^2 Xi exp 



-C(Xi) 



r sin(0) 



where Xi corresponds to the fraction of the photodissocia- 
tion rate due to radiation around the wavelength A^. For 
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HCN, 88% of the contribution comes from solar Lyma n a, 
i.e., i Ly a = 0.88 (jBockelee-Morvan fc Crovisierlll985h . To 
ease computations, we make the approximation 

This rough assumption is valid for small opacities (r > 
t thick) and if the water absorption cross-sections are the 
same at all wavelengths A^. Otherwise, it will slightly over- 
estimate the screening effect in the opaque region. We then 
define an effective water absorption cross-section for HCN 
(likewise for the other molecules) 

< cth 2 o^hcn >=^2xi(T K2 o{Xi) = 13.8 x 1(T 22 m 2 , (10) 

where most of the contribution co mes from Lyman a 
(a^dLy a) = 15 x 10" 22 m 2 , iLee fc Suto I . fl"98& 
iLewis et all . [1983). This approximation was validated by 
ourselves for HCN: we estimated numerically that using 
this mean value for the screening cross-section instead of 
summing over various wavelength intervals yields only a 
« 2% excess error in the estimate of the increase in the 
number of molecules due to screening. 

After solving the differential balance equation for the 
HCN density at the distance r from the nucleus in the coma, 
we find that 

, ,s <3hcn / r/3„ r \ 

n HC N(r, 0) = 2 exp -— G(- — ) (11) 

47rv cxp r y v CKp c HCN _ I _ w j 

with 

ChCN — < ^HaO^HCN > Qh 2 Q /-^ 
47TWexp 

The function G(x) = A f Q exp(— \/t)dt is equal to the ex- 
ponential integral G(x)= E2(l/a;), which can easily be com- 
puted by numerical integration. We note that G(oo) = 1.0, 
so that Eq. [TT] gives the classical Haser formula for negligi- 
ble opacities. 

We compared the production rates determined using the 
density from Eq. [TT] to those obtained with the Haser for- 
mula. The largest effect is for comet C/2006 PI with a « 
60% decrease of the HCN production rate, and ~ 40% de- 
crease for CH 3 OH, CH 3 CN, or HDO. We did not develop a 
full 3-D model to take into account phase angles different 
from 0° or 180°. But the comparison between the case of a 
phase angle of 180° , where we can simply use Eq. [TU and 
replacing Eq. U] by the averaged value in Eq. [5J only yields 
a 3% difference. 

Table [7] provides characteristic scalelengths for the three 
comets and for comparison C/1996 B2 (Hyakutake) and 
C/1995 01 (Hale-Bopp) observed in 1996-1997. In all cases, 
the optically thick region (at the Lyman a wavelength) is 
within the water coma (£h 2 o > (rthick{Ly a))) and if one 
takes into account OH, the assumption of an infinite life- 
time for the screening molecules is valid since the coma 
encompassing H2O and OH is definitely larger than the 
optically thick region and the HCN coma. One can also 
note that for comets C/1996 B2 and C/1995 01, L H cn > 
(rthick{Ly a)), so that the screening effect is negligible. 



4.3. Line shapes and gas expansion velocity 

We used the line sha pes to estimate the gas expansion ve- 
locity following e.g.. iBiver et all (| 19991 ). Table ED lists the 
velocities (VHM) corresponding to the half maximum in- 
tensity on both negative and positive sides of the lines. For 
each comet, we used the most reliably detected lines. To 
infer the gas velocity from the measured VHM , the pro- 
cesses contributing to line broadening should be considered. 
For the HCN(3-2) line, we took into account its hyper- 
fine structure (broadening of about 0.03 km s _1 and 0.06 
km s , on the negative and positive sides, respectively). 
Thermal broadening affects the widths of the lines, and de- 
pends on the actual gas temperature. When gas expansion 
velocities are in the range 1-2 km s , and temperatures 
are 100-300 K, thermal dispersion essentially smooths the 
line shape and widens it by less than 0.05 km s _1 . Optically 
thick lines, such as HCN(3-2) in C/2006 PI, are broadened 
by an additional 0.05 km s _1 when the expansion velocity 
is constant in the coma. 

At heliocentric distances < 0.25 AU, the photodis- 
sociation lifetimes of HCN, HNC, CH 3 OH, CS, HC 3 N, and 
H2CO are shorter than 1 h. In most cases, the beam size 
(5000 to 14000 km at the comet for C/2002 X5, C/2002 VI, 
and C/2006 PI near perihelion) is larger than the molecu- 
lar scale-length. Therefore, velocity measurements pertain 
to nucleus distances fixed by the molecular lifetimes and 
the line widths are mostly representative of the expansion 
velocity in the coma at a distance equal to the molecular 
scale-length. Figures I4T-HU1 show that the molecular lines 
have different widths. The measured "VHMs" velocities 
are plotted versus the molecular scale-lengths in Figs. [T3T 
1161 In the following sections, we constrain the evolution of 
the gas expansion velocity with distance to nucleus by fit- 
ting the measured VHMs, taking into account all sources 
of line broadening. 

4.3.1. Outgassing pattern and anisotropy 

Lines of C/2002 X5 detected around perihelion, and those 
of C/2006 PI, show an asymmetric shape, while C/2002 VI 
lines are relatively symmetric. 

For C/2002 X5, we likely observed a jet or outburst in 
progress on the night side. All the lines are strongly red- 
shifted (+0.25 km s" 1 on average, Table 0) suggesting gas 
outflow in a preferential direction, independent of accel- 
eration effects. The phase angle was small (26°, Table [T]), 
so the observed asymmetry corresponds to preferential out- 
gassing from the night side. Anisotropy in screening of pho- 
todissociation (Sect 4.2) also contributes to the asymmetry 
of the lines because lifetimes are longer on the night side. 
Considering the lines the less affected by screening (HC3N, 
CH3OH), we found that the outgassing rate is 1.35 ± 0.1 
times higher on the night side. The mean gas velocity is 
smaller (0.9 versus 1.2 km s" 1 ) on the dayside of the nu- 
cleus. This is most likely owing to shorter molecular life- 
times on the dayside, the measurements sampling thereby 
molecules closer to the nucleus (Sect. 4. 2). 

For C/2006 PI, the phase angle was large when the 
comet was observed (120-140°, Fig. [T2]). The line veloc- 
ity shifts are slightly negative (—0.05 km s _1 , Table 0]), 
which could indicate some excess emission from the night 
side. However, evidence of asymmetric gas acceleration is 
present. The mean gas velocities inferred from long-lived 
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Table 7. Typical scale lengths affected by screening 



Comet 


< r h > 


Qh 2 o 


^cxp 


Lh 2 o 


iOH 


< r t hick(Ly a) > 


£hcn 




[AU] 


[io 29 *- 1 ] 


[km s" 1 ] 


[km] 


[km] 


[km] 


[km] 


C/1996 B2 (Hyakutake) 


0.25 


5 


1.60 


7900 


14300 


90 


6600 


C/1995 Ol (Hale-Bopp) 


0.91 


100 


1.10 


70000 


128000 


2700 


58000 


C/2002 X5 (Kudo-Fujikawa) 


0.21 


35 


1.25 


3600 


7200 


820 


2900 


C/2002 VI (NEAT) 


0.12 


20 


2.00 


2000 


3800 


300 


1600 


C/2006 PI (McNaught) 


0.23 


300 


1.50 


8100 


14900 


5900 


4900 



Lh 2 o, Lou, Lhcn: unscreened scale lengths of H2O, OH, and HCN, taking into account solar activity, expansion velocity, and 
heliocentric distance. 

< rthickiJ^J Oi) >: mean radius of the optically thick Lyman a envelope (^C(Ly a)). 

Fig. 12. Sketch of the Earth-Comet-Sun 
plane, with the positions of each of the 
three comets at the epoch of the respec- 
tive observations. In addition (at a dif- 
ferent scale given by the 10000 km bar), 
three regions in the comae are plotted: 
in the shaded area, the region opaque to 
Lyman a (i.e. where opacity at 121.6 nm 
towards the Sun is larger than 1.0); in 
dashed line, the spherical region whose ra- 
dius corresponds to the Haser scale length 
of HCN (assuming a constant mean value 
for the expansion velocity) ; within the dot- 
ted lines, the region where less than 1- 
1/e of the HCN molecules are photodis- 
sociated (« equivalent to its scale length) 
when the screening is taken into account. 



10000 km 



% C/2006 PI (McNaught) 
\ k~ 16/01/2007 



L(HCN)= Vexp//S(HCN) - 
<L(HCN)> with screening 



Earth 



SUN 



1 • v-i ; C/2002 X5 (Kudo-Fujikawa) 
V.- 26/01/2003 



C/2002 VI (NEAT) 
17/02/2003 



-0.5 





Distance in AU 



0.5 



molecules (HCN, CH3OH) are 1.7 and 1.3 km s _1 on the 
rear (« Sun facing) and front side of the nucleus, respec- 
tively, but w 1.2 km s _1 on both sides for species with 
shorter lifetimes (CS and H 2 CO) (Table E). The less sig- 
nificant gas acceleration on the night side is likely due to a 
reduced photolytic heating caused by the screening of wa- 
ter photolysis. The small line blue-shifts suggest that the 
out gassing rate is barely larger (by a factor of 1.1 ± 0.1) 
in the day side, after taking into account screening effects 
(Sect. 4.2). 



In the case of C/2002 VI, the phase angle was close to 
90° (Table [lj, so that the line shapes are unaffected by 
day/night asymmetries. 



4.3.2. Expansion velocity acceleration 

Since there is evidence of gas acceleration in the coma, 
we used a variable expansion velocity in the coma 
Vex P (r), which is plotted in Figs. [TMTdI This function 
reproduces the general shape of gas-dynamic simulations 
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Fig. 13. The evolution of the expansion velocity with dis- 
tance to the nucleus in the coma of comet C/2002 X5 
(Kudo-Fujikawa) on 26.5 January 2003. The measured lines 
widths have been plotted with their errorbars at a dis- 
tance corresponding to the molecular lifetime. Arrows rep- 



Fig. 15. The evolution of the expansion velocity with dis- 
tance to the nucleus in the coma of comet C/2006 PI 
(McNaught) on 16.5 January 2007. Measurements of line 
widths from Table M are plotted as in Fig. [T3J The up- 



per values 



resent the effect on scalclcngth of the 0o = 1 3 x 10 J s larger velocities (v > 0). 
(Table 0) domain for CS. For HCN, HNC, and CH 3 OH, 
the arrow points to the increase in scalelength due to the 
screening. 

„ 2.5 



(r) and VHM) correspond to the ^dayside 



* HCN 

• CS (|8 =1-3.10 s ) 
■ CH,OH 



v elIp0 =lkm/s Xacc=B.l 




r h =0.135 AU 



10 10 10 10 

Distance from nucleus (molecules scale— length) (km) 



Fig. 14. The evolution of the expansion velocity with dis- 
tance to the nucleus in the coma of comet C/2002 VI 
(NEAT) on 16.5 February 2003. Measurements of line 
widths from Tableware plotted as in Fig. [TBI 



(jCombi and Smvthl . Il988f) . The parameter v CKp (r) Q is ad- 
justed to match the observed line widths (Table when 
computing line shapes. These velocity plots provide the two 
free parameters used to compute v oxp (r) and lead to the fol- 
lowing comments: 

1. C/2002 X5 at r h = 0.21 AU (Fig. [T3]): Positive and 
negative sides VHMs were both fitted using the same 
function. The computed line shapes for HC3N(28-27), 
CH 3 OH(5 - 5_i + 6 - 6_i + 7 - 7_i)E, and HCN(3- 
2), taking into account screening and a night/day out- 



xp (r) = v, 



cxpO 



0.6 + 0.3^/log(r/3) 



0.1 xacc (l + tanh (l.3Iog(=g^l))) 



+ VcxpO X 



, with r in km. 



a 



A HCN 

A HNC -, 
• CS QS =1-3.10 s ) 



v exp0 =0.9km/s Xacc=5.Z 




r h =0.S55 AU 



v expB =0.9km/s Xbcc=3.4 



10 10 10 10 

Distance from nucleus (molecules scale-length) (km) 



Fig. 16. The evolution of the expansion velocity with dis- 
tance to the nucleus in the coma of comet C/2006 PI 
(McNaught) on 17.5 January 2007. Scales as in Fig. [15] with 
upper values for dayside data. 



gassing rate ratio of 1.35, have exactly the measured 
width within the errorbars. Doppler shifts of the lines 
(+0.095, +0.172, +0.243 km s~\ predicted, respec- 
tively) are also in good agreement. 
2. C/2002 VI (Fig. [0}: Asymmetry in the lines (due to 
day /night out gassing asymmetry and screening effect) 
is neither seen nor expected because of the phase an- 
gle. We have thus averaged the VHM measurements 
on both sides of the lines . However, the two days of 
observations require two different velocity profiles, with 
stronger acceleration on the second day closer to the 
Sun. For this comet, we did not observe short-lived 
species (e.g., HC3N, H 2 CO, or H 2 S), so the velocity pro- 
file is poorly constrained close to the nucleus. The CS 
line widths suggest that the CS lifetime is significantly 
shorter than the HCN lifetime. 
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C/2006 PI: Data were divided into four subsets: posi- 
tive and negative velocity sides of the lines on the 16th 
(Fig.[H]) and 17th of January 2007 (Fig.[EHJ). The species 
with the smaller scale length (H2CO) was not observed 
on 17 Jan. but we assumed a similar acceleration in the 
coma for both days, scaled to the HCN lines VHMs. A 
stronger acceleration is present for molecules moving in 
an opposite direction to the observer (v > 0), which is 
consistent with stronger acceleration towards the day- 
side (Fig. [P2] phase angle of 120-140°< 90°). 



4.4. Constraints on CS and HNC lifetimes 



We assume that CS is the photodissociation product of CS2. 
It forms very close to the nucleus since CS2 photodissociates 
into CS in less than 30 seconds at rh < 0.25 AU. We assume 
that there is no addit i onal e xcitation effect to those not 
modeled in lBiver et al.l (|1999t ) that would strongly affect its 
rotational population mimicking photodissociation. Hence, 
the width of CS lines (clearly narrower than HCN lines; 
Figs @J HI En and Table O can be converted into a scale- 
length estimate using line-shape modeling with v CKp (r). 
The scale-length of CS is determined by the gas expan- 
sion velocity, its photodissociation rate, and the screen- 
ing effects. Since the photodissociation process of CS is 
unknown, we investigated /?o,cs = 1 _ 3 x 10 -5 s _1 for 
the photodissociation rate at 1 AU and cth 2 o^cs =0- 
12xl0~ 22 m 2 (i.e. values smaller than or comparable to 
HCN for the screening cross- section). Table M 

summarizes 

the findings. There is no unique solution, but grouping mea- 
surements comet by comet, the closer matches are found for 
(A>,cs > l-7x 10- 5 s-\ a H2 o^cs > 6x HT 22 m 2 ), (/3 ,cs = 
2.9 ± 0.7 x 10~ 5 s" 1 , for er H2 o^cs = 6 x 10~ 22 m 2 ), and 
(A>,cs = 4.3± 1.2 x 10- 5 s"\ (Th,o^cs > 6 x 10~ 22 m 2 ) for 
C/2002 X5, C/2002 VI, and C/2006 PI, respectively. The 
screening of photodissociation is necessary to achieve closer 
agreement between the night- and dayside measurements. 
In addition, the detection of CS at only 0.1 AU from the 
Sun in C/2002 VI requires a lifetime (/3 ,cs < 4x 10" 5 s" 1 ) 
long enough to get realistic abundances (CS/H 2 <C 1%, 
this ratio being < 0.2% at 1 AU in all comets where it has 
been measured). Values /?o,cs = 2.5 ± 0.5 x 10~ 5 s _1 and 
CTH2O-+CS = 6 x 10~ 22 m 2 are consistent with all measure- 
ments and are used to determine the production rates. 



Table 8. Constraints on the CS photodissociation rate at 
1 AU (A,cs) 



For HNC, we assume similar spectroscopic properties 
(e.g. UV absorption spectrum) to those of HCN. Slight 
differences in line shapes may underline differences. In 
C/2002 X5 and C/2006 PI, the HNC(3-2) line is slightly 
narrower than HCN(3-2) suggesting stronger photodisso- 
ciation, weaker screening, and/or other destruction pro- 
cesses in the coma. But this effect is not observed in 
C/2002 VI. A possible interpretation is that HNC is partly 
created in the coma by che mical reactions, as suggested by 
iRodeers fc Charnlevl (1998). This would partly compensate 
for the shortening of its scale-length. A parent scale-length 
of w 1300 km at rh = 0.2 AU, /?o,hnc ~ 1-3 x /?o,hcn, and 
(Th 2 0->hnc 1=3 0.9 x <7h 2 o->hcn would provide the closest 
agreements between model and data, but the constraints 
are not strong enough to infer a definite conclusion. 



< r h > Pocs 


CH 2 O^CS 


lifetime* 


[AU] [s- 1 ] 


[m 2 ] 


at 1AU [s] 


C/2002 X5: VHM > 


0.213 0.65 ± 0.15 x 10 _b 


no screening 


154000 


0.213 l.l±0.3xl0" 5 


3 x 10 22 


108100 


0.213 1.25 ± 0.35 x 10" 5 


6 x 10 22 


107500 


0.213 1.4±0.3xl0~ 5 


9 x 10 22 


106900 


0.213 1.6±0.3xl0" 5 


12 x 10 22 


104300 


C/2002 X5: VHM < 


0.213 2.7±0.7xl0 _b 


no screening 


37000 


0.213 2.3±0.7xl0~ 5 


3 x 10 22 


49400 


0.213 2.2±0.7xl0~ 5 


6 x 10 22 


55300 


0.213 2.1±0.7xl0~ 5 


9 x 10 22 


61000 


0.213 2.0±0.7xl0" 5 


12 x 10 22 


66600 


C/2002 VI: 16.6 Jan. 2003 


0.136 2.8±1.0xl0~ b 


no screening 


36000 


0.136 3.2±1.2xl0~ 5 


3 x 10 22 


40400 


0.136 3.8±1.6xl0" 5 


6 x 10 22 


40800 


C/2002 VI: 17.5 Jan. 2003 


0.108 2.05 ± 0.45 x 10 _b 


no screening 


49000 


0.108 2.45 ± 0.6 x 10~ 5 


3 x 10 22 


56400 


0.108 2.95 ± 0.75 x 10 -5 


6 x 10 22 


57700 


C/2006 PI: VHM < 


0.255 0.95 ± 0.4 x 10" b 


no screening 


105000 


0.255 2.05 ± 1.1 x 10~ 5 


3 x 10 22 


77000 


0.255 3.2±2.0xl0 -5 


6 x 10 22 


71000 


C/2006 PI: VHM > 


0.255 4.3±2.1xl0~ b 


no screening 


23000 


0.255 5.1 ±2.4xl0~ 5 


3 x 10 22 


36000 


0.255 6.4±3.4xl0 -5 


6 x 10 22 


34000 



* Determined from the ratio of the total number of molecules 
in the coma divided by the production rate, scaled as 1/r 2 to 1 
AU. 



5. Production rates and abundances 

The production rates were computed using models that 
incorporate collisions with neutrals and electrons, and ra- 
diative pum ping by the solar radiation ( Bi ver et al.l . [l999, 
I2000L l2006al ) . We did not consider infra-red pumping by the 
large and warm dust coma. Assuming a dust-to-gas ratio of 
1.0 (which might be underestimated for the dustier comets 
C/2006 PI and C/2002 VI), an d the simplified approach by 
ICrovisier and Encrenazl (|1983[ ). we estimate that for HCN 
or CH3OH in comet C/2006 PI (McNaught), the vibra- 
tional pumping by the dust infrared radiation is stronger 
than that by the solar radiation field within 2000-7000 km 
of the nucleus. However, collisions with water still dominate 
in a region twice as large and control the excitation of the 
rotational levels. 

We used the radial profiles of the gas expansion veloc- 
ity determined in the previous section. In some cases, we 
considered different velocity profiles to interpret the posi- 
tive and negative sides of the lines, but similar production 
rates were obtained as long as the velocity profiles provided 
a good fit to the line widths. To compute rotational level 
populations and collision rates with neutrals, we used the 
gas temperatures given in Section 4.1. Inferred production 
rates are given in Table [TO] As discussed in the next sec- 
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Table 9. Line widths and photodissociation rates 



UT date 


< r h > 


Line 




Beam 


size 


VHM 


V exp 


A) 




L d 


[mm/dd.dd] 


[AU] 






["] 


[km] 


[km s^ 1 ] 


[km s^ 1 ] 


[km s^ 1 ] 


Is" 1 ] 




[km] 


C/2002 X5: (2003) 


01/13.60 


0.553 


HCN(3-2) X 




9.6 


7190 


-1.09 ±0.07 


±1.20 ±0.07 


0.90 


1.7 x 10" 


"b 


14330 


01/26.54 


0.213 


HCN(3-2) 




9.4 


7990 


-0.97 ±0.03 


±1.42 ±0.02 


1.25 


1.92 x 10 


-5 


2950-4870 2 


01/26.54 


0.213 


HNC(3-2) 




9.2 


7820 


-0.86 ±0.11 


+1.28 ±0.03 


1.25 


1.92 x 10 


-5 


2950-4500 2 


01/26.56 


0.212 


CS(3-2) 




17.5 


14880 


-0.92 ±0.06 


+1.26 ±0.03 


1.20 


1-3x10" 


5 


1800-5400 


01/26.60 


0.212 


HC 3 N(28-27) 




9.7 


8250 


-0.73 ±0.18 


+0.93 + 0.11 


0.90 


6.8 x 10" 


■5 


690 2 


01/26.55 


0.213 


CH 3 OH((5 - 


5_ 


-i) + (6c 


,-6-1) 


+ (7 -7_i)E) 




















15.5 


13170 


-0.92 ±0.20 


+1.25 + 0.10 


1.25 


1.43 x 10 


-5 


4120-5130 2 


03/12.72 


1.184 


HCN(3-2) 




9.6 


7630 


-0.65 ±0.04 


+0.52 + 0.07 


0.65 


1.6 x 10" 


■5 


47460 


C/2002 VI: (2003) 


02/16.55 


0.136 


HCN(3-2) 




9.3 


6600 


-2.07 ±0.06 


±2.26 ± 0.04 


2.00 


1.8 x 10" 


-b 


2080-2950 2 


02/16.57 


0.135 


HNC(3-2) 




9.2 


6520 


-2.34 ±0.22 


+2.13 + 0.26 


2.00 


1.8 x 10" 


■5 


2030-2500 2 


02/16.55 


0.136 


CS(3-2) 




17.3 


12270 


-1.60 ±0.35 


+1.65 + 0.18 


1.60 


1-3x10" 


5 


1000-3000 


02/17.10 


0.120 


CH 3 OH((7 3 - 


7 2 )A+-+(6 3 - 6 2 )A+-+(5 3 - 5 2 )A+-+(4 3 - 4 2 )A+-+(7 3 


-7 2 )A"+) 














10.0 


7120 


-1.96 ±0.34 


+2.37 ± 0.43 


2.00 


1.38 x 10 


-5 


2090-2580 2 


02/17.50 


0.108 


HCN(3-2) 




9.3 


6650 


-2.07 ±0.07 


+2.24 + 0.06 


2.00 


1.8 x 10" 


■5 


1300-2250 2 


02/17.50 


0.108 


CS(3-2) 




17.3 


12370 


-2.60 ±0.60 


+1.65 + 0.10 


1.60 


1-3x10" 


5 


620-1870 


02/17.60 


0.105 


HCN(l-O) 




26.8 


19180 


-1.3 ±1.1 


+2.9 ± 1.1 


2.00 


1.8 x 10" 


■5 


1230-2200 2 


C/2006 PI: (2007) 


01/16.56 


0.228 


HCN(3-2) 




9.3 


5540 


-1.42 ±0.02 


+1.78 + 0.03 


1.50 


1.60 x 10 


-b 


4870-12300^ 


01/16.56 


0.228 


HCN(l-O) 




26.5 


15780 


-1.34 ±0.03 


+1.79 + 0.06 


1.50 


1.60 x 10 


—5 


4870-12300 2 


01/16.55 


0.228 


11200(312-2!! 


) 


10.8 


6430 


-1.10 ±0.20 


+1.10 + 0.15 


1.1 


2.0 x 10" 


-4 


290 3 -^1210 


01/16.56 


0.228 


CH 3 OH((4 - 


4_ 


-i) + (5c 


,-5-1) 


+ (6 - 6_i) ± (7 - 7_i)E) 


















15.5 


9230 


-1.45 ±0.15 


+1.90 + 0.30 


1.50 


1.32 x 10 


-5 


5900-10600 2 


01/17.56 


0.255 


HCN(3-2) 




9.3 


5620 


-1.33 ±0.05 


+1.67 + 0.05 


1.45 


1.60 x 10 


-5 


5890-12000 2 


01/17.55 


0.255 


HNC(3-2) 




9.2 


5480 


-1.23 ±0.11 


+1.60 + 0.15 


1.40 


1.60 x 10 


-5 


5670-10500 2 


01/17.56 


0.255 


CS(3-2) 




17.0 


10270 


-1.17 ±0.05 


+1.20 + 0.10 


1.2 


1-3x10" 


5 


2620-9800 


01/17.59 


0.256 


HCO+(1-0) 




26.2 


15830 


-2.80 ±0.53 


+2.03 + 0.53 











1 : Only low resolution (1 MHz = 1.2 km s" 1 ) spectra. 

2 : Scale- length taking into account screening (Sect. 4. 2): based on the equivalent lifetime that would lead to the same number of 
molecule in the coma without screening. 

''• When we assume that H2CO comes from an extended source with scale-length 1.75 x Z/h 2 co, this value increases by a factor 



4.23 (to get a 1/e decrease from the peak value at 1.30 x -Lh 2 co)- 



tion, they are possibly wrong by a factor 1.5, because of 
uncertainties in the model parameters. 



in these values by 50% affects the production rate determi- 
nations in the following way: 



5.1. Uncertainties due to model parameters 



The p hotodissociative lifetimes were taken from ICrovisierl 
(|1994|) . except for CS and HNC, which we constrained 
from the present observations (Sect. l4~4)l . We took into ac- 
count solar activity when computing the photodissociation 
rates (|Crovisierl I1989L 11994 iBockelee-Morvan fc Crovisie"rl 
Il985|) . The uncertainty introduced by our simplifying as- 
sumptions in the modeling of photolysis screening is not 
very large (< 15%, not considering water production rate 
uncertainties). Derived HCN production rates are a factor 
between 1.4 and 2.5 smaller than when this effect is omit- 
ted. 

The main source of error in the computation of produc- 
tion rates is the uncertainty in the gas kinetic temperature 
(except for C/2002 X5 at r h > 0.5 AU, where the tem- 
perature is well constrained). We assumed constant coma 
temperatures of 180, 150, and 300 K for near perihelion 
data of comets C/2002 X5, C/2002 VI, and C/2006 PI, 
respectively (Sect. 4.1). An increase (respectively decrease) 



Qhcn and Qhnc are increased by +33 to +40% 
(C/2006 PI) (respectively, decreased by 31% to 39%). 
the same trend is observed for Qcs and Qch 3 cn: ±45% 
for a ±50% variation in the gas temperature. 
Qch 3 oh is still more sensitive to the assumed tempera- 
ture, because of excitation of the torsional band at high 
T: a ±50% change in temperature causes changes of 
+51/-39%, +72/-53%, and +102/-64% in the produc- 
tion rates for C/2002 VI, C/2002 X5, and C/2006 PI, 
respectively. 

Qhc 3 n is minimal for T ss 200 K. The retrieved Qhc 3 n 
is 50% higher when either increasing or decreasing the 
nominal temperature (180 K) by 50%. 



Therefore, most abundances relative to HCN are not 
significantly sensitive to the assumed T. Two exceptions 
are the CH 3 OH/HCN (+20 to +40% for AT= +50%) and 
HC3N/HCN (+100% for AT = +50%) ratios. 
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1.5 10' 




26.4 



26.5 January 26.6 
Date (UT days) 



26.7 



Fig. 17. Molecular production rates in comet C/2002 X5 
(Kudo-Fujikawa) as a function of time on 26 January 2003: 
this plot shows the rapid increase over «3 h of observa- 
tions for HCN, CS, and HNC. The values for CH 3 OH and 
CH3CN are the averages over the full period. 



5.2. Short- and long-term variations in the production rates 

The most significant time variations are observed for 
C/2002 X5. Rapid variations are observed around perihe- 
lion (Fig. 1171) . Strong variations in th e light curve and w ater 
production rate are also reported bv lBout et al.l (|2003l) and 
iPovich et all (|2003l ). On 26.5 January 2003, the production 
rates of HCN, CS and HNC increased by ~ 70% over a time 
interval of « 0.11 days (Table [TOl Fig.[T7]). 

Over the 2-3 days of observations of comets C/2002 VI 
and C/2006 PI, production rates also varied. For 

C/2002 VI, an heliocentric dependence Qcs r h°' a to 

-4.5 



-1.5 



is observed, while for C/2006 PI Qhcn 



' h 

oc r. 



For C/2006 PI, the variation in the HCN production rate 
is steeper than the variation in the H2O production rate 



(2.4 x 10 



2!>, 



molec. s ) deduced from the OH observa- 



tions (Table|B]). This inconsistency could be due to rotation- 
induced time variations and the flux loss cause by anoma- 
lous refraction on 17 January being unable to be quantified. 

The long-term (over 2 months) evolution of the produc- 
tion rates of C/2002 X5 is shown in Fig. [H The HCN, 
HNC, and CH3OH production rates vary according to 
Q oc r - 3 - 5±0 - 5 ; p r e and post-perihelion, as for H 2 post- 
perihelion ((2.6 ±0.2) x io 2 8 r ^ 3 ' 5±0 ' 2 molec. s^ 1 ). 



5.3. Molecular abundances and discussion 

Table [TT] summarizes the molecular abundances relative to 
water and HCN in the three comets. 

The water production rates used as reference for the 
short heliocentric distances (r^ < 0.3 AU) are uncertain, 
as discussed in Section 3. For C/2002 X5, according to the 
short-term variations seen for all molecules, we assumed 

<3h 2 o = 35 + 10 x sin 



• fo t[d]-26.54\ 

m 12tt ' 2 J 



x 10 29 molec. s 1 . For 



C/2006 PI, we used Qh 2 o = 2.4 x 10 29 r^ J molec. s" 1 , 
and for C/2002 VI Qh 2 o = 2.7 x lO 29 ^ 1 molec. s" 1 , 
these values being possibly in error by ±50%. Resulting 
abundances relative to water are rough ly in agreement wit h 
abundances measured in other comets (jBiver et all l2002all . 
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Fig. 18. Molecular production rates in comet C/2002 X5 
(Kudo-Fujikawa) pre (left) and post (right) perihelion as 
a function of heliocentric distance. The water produc- 
tion r ates are based either on SOHO or O din observa- 
tions (jPovich et all 120031: iBiver et all l2007al) (black dots) 
or Nancay OH production rates multiplied by 1.1 (black 
squares). Dotted lines show the fitted evolutions (see text). 



The CH3OH/H2O ratio seems however quite low in comet 
C/2002 VI and especially in C/2006 PI. The HDO/H 2 
ratio in comet C/2006 PI, is not well constrained, both 
OH and HDO being only marginally detected, but produc- 
tion rate s are compatible with t h e rati o measured in other 
comets (iBockelee-Morvan et all . Il998t iMeier et all 119981: 
Uehin et all I2009L » 6 x 1Q- 4 1. 

We now turn our discussion to abundances relative to 
HCN, which are more reliable since HCN was always ob- 
served simultaneously. The main observed features are: 

— A decrease in the CH3OH/HCN ratio at low r/j oc r" 5 is 
observed for C/2002 X5, confirming the trend observe d 
in comet C/1996 B2 (Hyakutake) (|Biver et all 11999(1 . 
Using this 5 variation, we infer CH3OH/HCN ratios 
at 1 AU of 18 and 10 for C/2002 VI and C/2006 PI, 
respectively, which are within the range of abundances 
measured in comets, tho ugh C/2006 PI belo ngs to the 
methanol-poor category (jBiver et all l2002aT) 



-1.0 



A simi- 



The CS/HCN ratio varies as 0.35 x 
lar heliocentric dependence (^°' 8 ) was measured in 
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Table 10. Molecular production rates 



UT date Production rate in lO^molec. s 



[mm/dd.d] [AU] HCN HNC CH 3 CN CH 3 OH HaCO 1 CO CS HC 3 N PCS 

C/2002 X5 (Kudo- Fujikawa) (2003) 



01/13.61 


0.553 


1.6 ±0.3 


0.32 ±0.10 




25 ±8 




< 95 2 








01/26.47 


0.214 


56.6 ±3.0 


11.3 ±2.4 










101.2 ± 8.7 






01/26.51 


0.214 


63.8 ±3.6 


14.5 ± 1.4 










80.1 ±8.3 






01/26.54 


0.213 


88.7 ±5.2 


19.2 ±0.8 


5.6 ± 2.2 


841 ± 82 






130.8 ± 5.3 






01/26.58 


0.212 


95.5 ±1.6 


16.1 ± 1.2 
















01/26.61 


0.212 














156.0 ±6.1 


34.8 ± 5.6 




03/12.72 


1.184 


0.19 ±0.02 


< 0.064 


< 0.15 


5.0 ± 1.3 






< 0.19 






C/2002 VI (NEAT) (2003) 


02/16.6 


0.135 


34.0 ± 1.0 


10.0 ±0.9 


< 5.6 


220 ± 60 






102 ± 8 




< 460 


02/17.5 


0.108 


41.1 ±2.5 






261 ± 50 


< 389 


< 1075 


144 ± 27 


< 10.3 


< 581 


C/2006 PI (McNaught) (2007) 
HCN 


HNC 


CH 3 CN 


CH3OH 


H2CO 1 


CO 


CS 


HCOOH 


HDO 


01/15.65 


0.207 


359 ± 42 










< 37500 








01/16.54 


0.228 


286 ± 28 






1450 


590 






222 ± 60 3 


90 


01/16.58 


0.229 


242 ± 23 






±100 


±110 








±30 3 


01/17.55 


0.255 


146 ± 14 


16.7 ±2.7 


10.6 ±4 


645 ± 77 






217 ±33 






01/17.58 


0.256 


142 ± 15 










4290 3 









±1380 



: Assuming that H2CO comes from an extended source with scale-length 1.75 x Lh 2 co- 
2 : Uncertainty in the tuning: the backend may not have been properly connected to the receiver. 
3 : Marginal values. 



comet C/1996 B2 dBiver et all Il999t). In comet 153P, 
CS/HCN followed 0.55r^ 07 (|Biver et all . l2006ajh The 
puzzling heliocentric variation in the CS/HCN ratio is 
confirmed down to very low rv 
— Measuring the HNC/HCN ratio at small was one 
of the main objectives of our analysis of these obser- 
vations. HNC was easily detected in the three comets 
with an abundance ratio HNC/HCN between 0.11 and 
0.29 at r/j < 0.26 AU. The measurements in comets 
C/2002 X5 and C/2002 VI were included in a previous 
study of the heliocentric variation of HNC/HCN based 
on a sample of 11 c omets at in the range 0.14-1.5 
AU (|Lis et all . l2008li . The HNC/HCN sample was fit- 
ted by a power law in r^ 2 ' 3 , with an indication being 
found of a possible flattening at < 0.5 AU. The values 
inferred for our three comets do not show a trend for 
the more productive comets being enriched in HNC, as 
expected for a format i on of HNC by chemical reactions 
(Rod gers fc Charnlevl |1998[ ). but rather the opposite. 
The HNC/HCN ratio in C/2006 PI is 0.13, smaller than 
the value in comet Hale-Bopp at 1 AU (« 0.25), this 
comet being slightly more active than comet Hale-Bopp. 
Possible ori gins of HNC in cometary at mosp heres are 
discus sed in iRodgers &; Charnlevl (|l998h and iLis et all 
( 2008). Destruction of HNC by reaction processes pos- 
sibly takes place at low r^. Small differences in the line 
shapes between HCN and HNC suggest that HNC might 
be partly produced in the coma and destroyed further 
away from the nucleus or be more sensitive to photodis- 
sociation than HCN. 

Comet C/2006 PI has compositional similarities with 
fragments B and C of comet 73P/Schwassmann-Wachmann 
3. All are CH 3 OH and CO-poor, a nd the H 2 C O /HCN 
CH3CN/HCN ratios are similar (jBiver et all . l2006ri 



iDello Russo et all . l2007bf ). However, C/2006 PI is rich in 
volatile hydrocarbons and in NH 3 , while 73P is d e pleted in 
these compounds according to lDello Russo et"al~l ( 2007aD . 

Comet C/2002 X5, on the other hand, has a normal 
CH3OH abundance. It is unusually rich in HC3N. The 
HC3N /HCN ratio is four time s high er than in comet Hale- 
Bopp (jBockelee-Morvan et al.l . l2000D . and still higher tha n 
upper limits found in some comets (jBiver et all [2006a). 
Interest ingly, C/2002 X5 be longs to the class of carbon-rich 
comets (jPovich et all 120031 ) . Finally, the large abundances 
of HC 3 N, HNC, and CH 3 CN relative to HCN in this comet 
imply that these species are significant contributors to the 
production of CN radicals in this comet. 

5.4. Search for refractory molecules 

SiO was searched through its J(6-5) transition at 
260518.020 MHz in comets C/2002 X5 and C/2002 VI 
around perihelion (Tables I2"l3"j) . It was not detected down 
to about 0.01% relative to water, assuming thermal equi- 
librium and a photodissociation rate of 10~ 5 s . The same 
tuning at 260 GHz also allowed us to search for NaCl J(20- 
19) line in C/2002 VI, while the KC1 J(30-29) line was 
searched for at 230 GHz in C/2002 VI and C/2006 PI, in 
both case without success. 



6. Conclusion 

We have reported on our analysis of millimeter spectro- 
scopic observations of three bright comets (C/2002 X5 
(Kudo-Fujikawa), C/2002 VI (NEAT), and C/2006 PI 
(McNaught)) when they were close to the Sun (r^ < 
0.25 AU). Our main challenge has been to cope with 
ephemeris uncertainties by searching for the peak of bright- 
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Table 11. Relative abundances 



Molecule 


C/2002 X5 (Kudo-Fujikawa) 


C/2002 VI (NEAT) 


C/2006 PI (McNaught) 


< r h > i 


0.5-1.2 


0.21 


0.12 


0.24 


Abundance relative to water 


HCN 


0.11 ±0.02% 


0.23 ± 0.02% 


0.17 ±0.01 % 


0.13 ±0.02 % 


HNC 


0.027-< 0.03% 


0.047 ± 0.008% 


0.05 ± 0.005% 


0.011 ±0.002% 


CH 3 CN 




0.016 ± 0.006% 


< 0.024 % 


0.008 ± 0.003% 


HC 3 N 




0.081 ± 0.013% 


< 0.045 % 




CH3OH 


2.3 ± 0.5 % 


2.4 ± 0.2% 


1.1 ±0.2 % 


0.6 ± 0.2% 


H 2 CO 






< 1.95 % 


0.29 ± 0.06% 


CO 


< 8 % 




< 4.3 % 


3± 1% 


HCOOH 




< 0.9 % 




0.11 ±0.03% 


HDO 








0.045 ± 0.015% 


cs 


< 0.088% 


0.36 ± 0.02% 


0.54 ± 0.08% 


0.15 ±0.02% 


ocs 






< 1.54 % 




so 






< 0.81 % 




SiO 




< 0.008 % 


< 0.011% 




Abundance relative to HCN 


HNC 


0.20 ± 0.07-< 0.33 


0.198 ±0.024 


0.29 ±0.028 


0.114 ±0.026 


CH 3 CN 




0.07 ±0.03 


< 0.14 


0.08 ±0.03 


HC 3 N 




0.36 ± 0.06 


< 0.25 




CH3OH 


20 ±5 


11 ± 1 


6.4 ± 1.0 


5.1 ±0.3 


H 2 CO 






< 9.5 


2.2 ± 0.4 


CO 


< 60 




< 26 


31 ± 10 


HCOOH 




< 4 




0.8 ±0.2 


HDO 








0.3 ±0.1 


CS 


< 0.96 


1.55 ±0.17 


3.2 ±0.3 


1.51 ±0.24 


OCS 






< 8.6 




SiO 




< 0.04 


< 0.068 





ness in HCN emission using coarse mapping. We have ob- 
tained the following results: 

— The screening effect in photolytic processes, which af- 
fects the distribution of the molecules in the coma, was 
modeled to determine accurate production rates. This 
effect is significant, especially for C/2006 PI where it 
increases by 40 to 60% the number of molecules in the 
coma. 

— The various lines have different widths. Velocity varia- 
tions in the coma were considered to interpret the line 
shapes. This provided constraints on the CS and HNC 
scale-lengths. 

— The CS photodissociation rate was estimated to 2.5 ± 
0.5 x 10 -5 s _1 ( rh = 1 AU), which i s compatible with 
other estimates dBoissier et all [20071 ). 

— The HNC photodissociation lifetime is found to be 
slightly shorter (by « 30%) than the lifetime of HCN, 
though this estimate does not consider destruction 
paths of HNC by chemical reactions. 

— The CS/HCN production rate ratio in cometary atmo- 
spheres follows a heliocentric dependence ~ r h°' S down 
to r/j = 0.1 AU. We rule out artefacts related to un- 
certainties in CS lifetime. The origin of this variation 
remains unexplained. 

— The CH3OH/HCN production rate ratio decreases with 
decreasing according to r°' 5 . 

— HNC/HCN ratios are high (0.11 to 0.29) in the range 
r h = 0.14-0.25 AU, consistent with HNC being a by- 
product of the thermal degradation of organic grains 
(|Lis et all . I2008D . A lower HNC abundance is measured 
for the more productive comet, possibly indicating de- 



struction of this reactive molecule by chemical pro- 
cesses. 

— The H2CO/HCN abundance ratio measured in 
C/2006 PI (2.2 at 0.23 AU) and C/2002 VI (< 9.5 at 
0.12 AU) is in the ra nge of values measured in c omets at 
1 AU from the Sun (iBiver et all . l2002ai l2006al 1.6-10). 
This possibly rules out H2CO production fro m the ther- 
mal d egradation of polymers as proposed bv lFrav et al.l 
(120061) . 

— The HC 3 N abundance in C/2002 X5 is higher, by 
a factor of four or more, than in any of the other 
eight comets in which it had previously been measured. 
HC 3 N, CH 3 CN, and HNC are all significant contribu- 
tors to the production of CN radicals in this comet. 

— Searches for SiO, NaCl, and KC1 were unsuccessful. The 
upper limit set for SiO is < 10 -4 relative to H2O. 
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Table 2. Molecular observations in comet C/2002 X5 (Kudo-Fujikawa) 



UT date (2003) 


< r h > 


< A > 


Int. time 


Line 




J' T b dv 


Velocity offset 


offset 


[mm/dd.dd-dd.dd] 


[AU] 


[AU] 


[min] 






[K km s" 1 ] 


[km s^ 1 ] 




IRAM 30-m: 


01/13.55-13.66 


0.553 


1.032 


60 


HCN(l-O) 




0.119 ±0.024 


-0.15 ±0.25 


3.5" 


01/13.55-13.66 


0.553 


1.032 


40 


HCN(3-2) 




1.405 ±0.058 


±0.02 ± 0.08 


3.5" 


01/26.46-26.46 


0.214 


1.172 


10 


HCN(3-2) 




1.075 ±0.190 


±0.22 ± 0.22 


10" 


01/26.47-26.49 


0.214 


1.172 


10 


HCN(3-2) 




4.640 ± 0.190 


±0.24 ±0.05 


6" 


01/26.50-26.50 


0.214 


1.172 


5 


HCN(3-2) 




6.046 ± 0.266 


±0.29 ± 0.09 


4.5" 


01/26.51-26.57 


0.213 


1.172 


90 


HCN(3-2) 




14.268 ± 0.084 


±0.24 ±0.01 


1.0" 


01/26.58-26.58 


0.213 


1.172 


5 


HCN(3-2) 




11.940 ±0.335 


±0.30 ± 0.04 


4.0" 


03/12.65-12.79 


1.184 


1.096 


90 


HCN(3-2) 




0.275 ±0.027 


-0.16 ±0.04 


3.0" 


01/13.55-13.66 


0.553 


1.032 


20 


HNC(3-2) 




0.326 ±0.101 


±0.36 ± 0.40 


3.5" 


01/26.46-26.46 


0.214 


1.172 


10 


HNC(3-2) 




< 0.501 




13" 


01/26.47-26.49 


0.214 


1.172 


10 


HNC(3-2) 




0.660 ± 0.140 


±0.40 ± 0.26 


6.8" 


01/26.50-26.57 


0.213 


1.172 


90 


HNC(3-2) 




2.415 ± 0.093 


±0.33 ±0.15 


2.8" 


01/26.58-26.58 


0.213 


1.172 


5 


HNC(3-2) 




1.494 ±0.399 


±0.02 ± 0.39 


4.4" 


03/12.65-12.79 


1.184 


1.096 


90 


HNC(3-2) 




< 0.106 




3.0" 


01/26.47-26.63 


0.213 


1.172 


135 


CH 3 CN(8,0-7,0) 




0.068 ± 0.016 


±0.29 ±0.25 


1.6" 










CH 3 CN(8,1-7,1) 




0.041 ±0.015 


-0.08 ± 0.33 












CH 3 CN(8,2-7,2) 




0.028 ± 0.019 














CH 3 CN(8,3-7,3) 




0.091 ±0.019 


-0.03 ±0.27 




03/12.65-12.79 


1.184 


1.096 


90 


CHaCN^) 1 




< 0.054 




3.0" 


01/26.60-26.63 


0.212 


1.172 


25 


HC 3 N(28-27) 




0.700 ±0.112 


-0.03 ±0.16 


1.5" 


01/13.55-13.66 


0.553 


1.032 


60 


CO(2-l) 




< 0.094 1 




3.5" 


01/26.46-26.46 


0.214 


1.172 


10 


CS(3-2) 




0.278 ± 0.052 


±0.45 ± 0.24 


10" 


01/26.47-26.49 


0.214 


1.172 


10 


CS(3-2) 




0.407 ± 0.057 


±0.23 ±0.18 


5.5" 


01/26.51-26.51 


0.214 


1.172 


5 


CS(3-2) 




0.463 ± 0.048 


±0.04 ±0.13 


1.5" 


01/26.52-26.57 


0.213 


1.172 


85 


CS(3-2) 




0.744 ± 0.030 


±0.31 ±0.05 


1.3" 


01/26.50-26.58 


0.213 


1.172 


10 


CS(3-2) 




0.787 ± 0.058 


±0.10 ±0.09 


3.7" 


01/26.60-26.63 


0.212 


1.172 


25 


CS(3-2) 




0.875 ± 0.034 


±0.25 ±0.05 


1.5" 


03/12.65-12.79 


1.184 


1.096 


135 


CS(3-2) 




< 0.032 




3.0" 


01/26.60-26.63 


0.212 


1.172 


25 


SiO(6-5) 




< 0.516 




1.5" 


01/13.55-13.66 


0.553 


1.032 


60 


CH 3 OH(1,0-1,-1 


E 


0.026 ±0.013 


-0.51 ± 0.36 


3.5" 








CH 3 OH(2,0-2,-l 


E 


0.065 ± 0.018 


±0.58 ± 0.39 












CH 3 OH(3,0-3,-l 


E 


0.099 ±0.018 


-0.33 ±0.21 












CH 3 OH(4,0-4,-l 


E 


0.049 ± 0.018 


-0.25 ±0.54 












CH 3 OH(5,0-5,-l 


E 


0.065 ± 0.018 


±0.03 ± 0.34 












CH 3 OH(6,0-6,-l 


E 


0.047 ± 0.018 


±0.25 ± 0.38 












CH 3 OH(7,0-7,-l 


E 


0.031 ±0.018 


±0.46 ± 0.64 




01/26.47-26.63 


0.213 


1.172 


135 


CH 3 OH(1,0-1,-1 


E 


0.072 ± 0.027 


-0.16 ±0.38 


3.0" 








CH 3 OH(2,0-2,-l 


E 


0.101 ±0.030 


±0.75 ± 0.41 












CH 3 OH(3,0-3,-l 


E 


0.142 ± 0.028 


-0.28 ± 0.22 












CH 3 OH(4,0-4,-l 


E 


0.113 ±0.030 


-0.11 ± 0.31 












CH 3 OH(5,0-5,-l 


E 


0.197 ± 0.028 


±0.30 ±0.17 












CH 3 OH(6,0-6,-l 


E 


0.197 ± 0.029 


±0.13 ±0.17 












CH 3 OH(7,0-7,-l 


E 


0.214 ± 0.028 


-0.10 ±0.15 




03/12.65-12.79 


1.184 


1.096 


135 


CH 3 OH(1,0-1,-1 


E 


0.042 ± 0.012 


±0.24 ±0.12 


3.0" 










CH 3 OH(2,0-2,-l 


E 


0.024 ± 0.014 














CH 3 OH(3,0-3,-l 


E 


0.011 ±0.013 














CH 3 OH(4,0-4,-l 


E 


0.015 ±0.014 














CH 3 OH(5,0-5,-l 


E 


0.043 ± 0.016 


-0.12 ± 0.15 












CH 3 OH(6,0-6,-l 


E 


0.036 ±0.012 


-0.03 ±0.12 












CH 3 OH(7,0-7,-l 


E 


0.027 ±0.012 


-0.03 ±0.16 












sum of 7 lines 




0.150 ±0.037 


-0.01 ±0.09 





1 Sum of CH 3 CN(8,0-7,0), CH 3 CN(8,1-7,1) and CH 3 CN(8,2-7,2) lines. 
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Table 3. Molecular observations in comet C/2002 VI (NEAT) 



UT date (2003) 


< r h > 


< A > 


Int. time 


Line 


J T b dv 


Velocity offset 


offset 


[mm/dd.dd-dd.dd] 


[AU] 


[AU] 


[min] 




[K km s^ 1 ] 


[km s^ 1 ] 




IRAM 30-m: 


02/16.49-16.59 


0.136 


0.978 


50 


HCN(3-2) 


4.146 ±0.061 


-0.02 ± 0.03 


1.8" 




0.136 


0.978 


19 


HCN(3-2) 


2.173 ±0.112 


-0.01 ± 0.11 


5.1" 


02/17.48-17.52 


0.108 


0.985 


19 


HCN(3-2) 


3.558 ± 0.069 


±0.01 ± 0.04 


1.8" 


02/17.48-17.50 


0.109 


0.985 


9 


HCN(3-2) 


1.819 ±0.140 


-0.18 ±0.17 


5.7" 


02/17.55-17.66 


0.106 


0.987 


85 


HCN(l-O) 


0.116 ±0.021 


±0.06 ± 0.32 


1.8" 


02/16.49-16.59 


0.135 


0.978 


29 


HNC(3-2) 


0.648 ±0.113 


±0.09 ± 0.35 


5.5" 


02/16.49-16.66 


0.134 


0.978 


85 


HNC(3-2) 


0.914 ± 0.060 


±0.25 ±0.13 


2.8" 


02/16.49-16.59 


0.135 


0.978 


27 


CS(3-2) 


0.242 ± 0.054 


±0.57 ± 0.42 


5.5" 


02/16.49-16.66 


0.134 


0.978 


87 


CS(3-2) 


0.336 ± 0.026 


±0.05 ±0.13 


2.8" 


02/17.48-17.52 


0.108 


0.985 


22 


CS(3-2) 


0.373 ± 0.046 


-0.34 ±0.21 


1.7" 


02/17.48-17.50 


0.109 


0.985 


6 


CS(3-2) 


0.180 ±0.085 




5.6" 


02/16.49-17.50 


0.130 


0.980 


22 


CH 3 CN(8-7) 1 


< 0.126 




2.5" 


02/17.61-17.66 


0.106 


0.987 


35 


CO(2-l) 


< 0.161 




2.8" 


02/16.54-17.66 


0.120 


0.982 


80 


UCb(li-ll) 


< 0.029 




2.1" 








HC 3 N(16-15) 
H 2 CO(2 2 -loi) 


< 0.031 
0.032 ±0.012 




2.1" 
2.1" 


02/17.55-17.60 


0.107 


0.986 


50 


H 2 CO(3i2 - 2n) 


< 0.136 




1.0" 


02/16.60-17.60 


0.120 


0.982 


107 


CH 3 OH(3,3-3,2)A ± 2 

CH 3 OH(4,3-4,2)A ± 

0113011(5,3-5,2^ 

CHaOH^^-e^^A* 

CHgOH^^-^A* 

0113011(8,3-8,2^ 


0.072 ± 0.045 
0.213 ±0.046 
0.092 ± 0.044 
0.182 ±0.047 
0.198 ±0.041 
0.057 ±0.041 


±0.20 ± 0.35 

±0.31 ± 0.48 
-0.28 ± 0.29 


1.8" 


02/16.60-17.60 


0.120 


0.982 


107 


SO(5 6 -4 s ) 


< 0.071 




1.8" 


02/17.61-17.66 


0.106 


0.987 


35 


SiO(6-5) 


< 0.190 




1.7" 



1 : sum of CH 3 CN(8,0-7,0), CH 3 CN(8,1-7,1), CH 3 CN(8,2-7,2) and CH 3 CN(8,3-7,3) lines.; 
2 : sums of twin lines; 
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Table 4. Molecular observations in comet C/2006 PI (Me Naught) 



UT date (2007) 


< r h > 


< A > 


Int. time 


Line 


J T b dv 


Velocity offset 


offset 


[mm /dd . dd-dd .dd] 


[AU] 


[AU] 


[min] 




[K km s^ 1 ] 


[km s^ 1 ] 




IRAM 30-m: 


01/15.65 


0.207 


0.817 


0.1 


HCN(3-2) 


83.74 ± 7.63 


-0.27 ±0.16 


2.6" 


0.207 


0.817 


0.2 


HCN(3-2) 


57.66 ± 5.82 


+0.33 ±0.18 


5.0" 




0.207 


0.817 


2.0 


HCN(3-2) 


31.78 ±3.07 


±0.40 ±0.18 


7.4" 




0.207 


0.817 


0.3 


HON (3-2) 


22.89 ±5.08 


-0.21 ±0.38 


8.9" 




0.207 


0.817 


0.8 


HCN(3-2) 


13.76 ±3.32 


±0.33 ± 0.44 


12.1" 


01/16.53 


0.228 


0.821 


8.0 


HCN(3-2) 


48.45 ± 0.42 


-0.07 ±0.02 


4.9" 


01/16.54 


0.228 


0.821 


8.0 


HCN(3-2) 


63.65 ± 0.39 


-0.09 ±0.01 


2.1" 


01/16.55 


0.228 


0.821 


20.0 


HCN(3-2) 


37.47 ±0.26 


-0.06 ±0.02 


6.5" 


01/16.54 


0.228 


0.821 


4.0 


HON (3-2) 


27.28 ±0.57 


-0.03 ±0.05 


8.2" 


01/16.58 


0.229 


0.822 


12.0 


HCN(3-2) 


47.89 ± 0.28 


-0.03 ±0.01 


4.7" 


01/16.59 


0.229 


0.822 


8.0 


HON (3-2) 


16.60 ±0.33 


±0.04 ± 0.05 


9.3" 


01/17.55 


0.255 


0.832 


2.0 


HCN(3-2) 


35.91 ±0.59 


-0.05 ±0.04 


3.3" 


01/17.54 


0.255 


0.832 


10.0 


HCN(3-2) 


21.81 ±0.28 


-0.02 ±0.03 


5.3" 


01/17.56 


0.255 


0.833 


31.0 


HCN(3-2) 


15.58 ±0.16 


±0.04 ± 0.02 


7.1" 


01/17.55 


0.255 


0.832 


5.0 


HON (3-2) 


14.44 ± 0.40 


-0.02 ±0.06 


8.5" 


01/17.58 


0.256 


0.833 


20.0 


HON (3-2) 


6.72 ±0.22 


-0.04 ±0.07 


10.7" 


01/17.58 


0.256 


0.833 


12.0 


HCN(3-2) 


5.69 ±0.27 


-0.06 ±0.11 


12.4" 


01/16.54 


0.228 


0.821 


20.0 


HCN(l-O) 


1.723 ±0.057 


±0.16 ±0.05 


3.9" 


01/16.55 


0.228 


0.821 


16.0 


HCN(l-O) 


1.505 ± 0.066 


±0.02 ±0.06 


6.8" 


01/16.54 


0.228 


0.821 


4.0 


HCN(l-O) 


1.574 ±0.131 


-0.13 ±0.11 


8.1" 


01/16.59 


0.229 


0.822 


12.0 


HCN(l-O) 


1.920 ±0.031 


-0.01 ±0.02 


4.8" 


01/16.59 


0.229 


0.822 


8.0 


HCN(l-O) 


1.261 ±0.037 


-0.08 ±0.03 


10.2" 


01/17.55 


0.255 


0.832 


2.0 


HNC(3-2) 


4.05 ±0.55 


-0.06 ±0.18 


3.7" 


01/17.54 


0.255 


0.832 


9.0 


HNC(3-2) 


3.16 ±0.24 


±0.03 ± 0.09 


4.7" 


01/17.56 


0.255 


0.833 


23.0 


HNC(3-2) 


1.53 ±0.16 


±0.11 ±0.12 


6.7" 


01/17.54 


0.255 


0.832 


14.0 


HNC(3-2) 


1.52 ±0.20 


±0.13 ±0.16 


9.1" 


01/17.54 


0.255 


0.832 


8.0 


HNC(3-2) 


0.52 ±0.26 




11.5" 


01/17.55 


0.255 


0.832 


2.0 


CS(3-2) 


2.84 ±0.20 


±0.07 ±0.08 


3.4" 


01/17.54 


0.255 


0.832 


8.0 


CS(3-2) 


2.47 ±0.10 


-0.04 ±0.04 


4.7" 


01/17.56 


0.255 


0.833 


24.0 


CS(3-2) 


1.53 ±0.06 


-0.03 ±0.04 


6.8" 


01/17.54 


0.255 


0.832 


13.0 


CS(3-2) 


1.17 ±0.08 


-0.09 ±0.08 


9.1" 


01/17.54 


0.255 


0.832 


9.0 


CS(3-2) 


1.01 ±0.09 


±0.11 ±0.10 


11.4" 


01/17.51-17.58 


0.255 


0.832 


47 


CH 3 CN(8,0-7,0) 


0.148 ± 0.040 


±0.45 ±0.20 


6.9" 










CH 3 CN(8,1-7,1) 


0.080 ± 0.044 


-0.60 ±0.47 












CH 3 CN(8,2-7,2) 


0.004 ± 0.039 














CH 3 CN(8,3-7,3) 


0.104 ± 0.042 


-1.01 ±0.62 












Sum of the 4 lines 


0.334 ± 0.083 






01/16.52-16.60 


0.228 


0.821 


48 


CH 3 OH(1,0-1,-1)E 


0.203 ± 0.050 




4.9" 








CH 3 OH(2,0-2,-l)E 


0.160 ±0.047 


-0.52 ±0.44 












CH 3 OH(3,0-3,-l)E 


0.309 ± 0.050 














CH 3 OH(4,0-4,-l)E 


0.327 ± 0.048 


-0.06 ±0.21 












CH 3 OH(5,0-5,-l)E 


0.305 ± 0.047 


±0.14 ±0.22 












CH 3 OH(6,0-6,-l)E 


0.347 ± 0.047 


-0.05 ±0.20 












CH 3 OH(7,0-7,-l)E 


0.297 ± 0.042 


-0.24 ±0.20 




01/16.54-16.60 


0.228 


0.821 


8 


CH 3 OH(J,0- J.-lJE 1 


0.857 ±0.198 


-0.60 ±0.35 


8.9" 


01/17.51-17.58 


0.255 


0.832 


48 


CH 3 OH(1,0-1,-1)E 


0.024 ± 0.035 




6.7" 










CH 3 OH(2,0-2,-l)E 


-0.021 ±0.039 














CH 3 OH(3,0-3,-l)E 


0.081 ± 0.035 














CH 3 OH(4,0-4,-l)E 


0.028 ± 0.042 


±0.64 ±0.60 












CH 3 OH(5,0-5,-l)E 


0.095 ± 0.040 


-0.47 ±0.46 












CH 3 OH(6,0-6,-l)E 


0.103 ±0.045 


-0.17 ±0.47 












CH 3 OH(7,0-7,-l)E 


0.082 ± 0.040 


-0.55 ±0.40 












CH 3 OH(J,0- J.-lJE 1 


0.380 ± 0.079 


-0.40 ±0.29 


6.7" 


01/16.52-16.60 


0.228 


0.821 


60 


H 2 CO(3i2 - 2u) 


0.334 ± 0.050 


-0.18 ±0.20 


6.1" 


01/16.52-16.59 


0.228 


0.821 


36 


HDO(3i2 - 2u) 


0.227 ± 0.064 


-0.68 ±0.40 


4.2" 


01/15.60 


0.207 


0.817 


5.0 


CO(2-l) 


< 2.31 




12.1" 


01/17.59-17.61 


0.256 


0.833 


24 


CO(2-l) 


0.30 ± 0.09 


-0.07 ±0.43 


11.2" 


01/17.59-17.61 


0.256 


0.833 


24 


HCO+(1-0) 


0.36 ±0.05 


-0.05 ±0.27 


10.9" 



1 : sum of CH 3 OH lines J, K = (4,0-4,-l)E, (5,0-5,-l)E, (6,0-6,-l)E and (7,0-7,-l)E. 



